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ABSTRACT 

The twist-four quark and gluon distributions in a transversely polarized nucleon are 
identified and relations among them are discussed by using QCD equations of motion. Their 
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nucleon-nucleon collisions is calculated. 
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I. INTRODUCTION 

There has been considerable interest recently in the polarization-dependent quark distri- 
butions of a transversely polarized nucleon.^"^ These distributions, which are on the same 
footing as those of an unpolarized or longitudinally polarized nucleon, are just as important 
for characterizing the nucleon's high-energy structure. Being universal quantities — meaning 
that they describe the nucleon rather than a particular process — these distributions can 
be measured in a variety of experiments involving large momentum transfers. Although at 
present no parton distribution can be calculated ab-initio from QCD, nonetheless this may 
eventually become possible. 

A very interesting distribution, hi{x), called the transversity distribution in ref. 2, was 
first identified by Ralston and Soper. That this involves a helicity flip of the participating 
quark was recognized by Artru and Mekhfi.-*^ The hi{x) is a twist- two quantity, meaning that 
it enters into the cross section of a hard process, such as Drell-Yan lepton pair production, 
unsuppressed by inverse powers of the hard momentum Q . It has been proposed that hi (x) 
be measured using transversely polarized beams at RHIC.^ The O o^" twist-four, cor- 

rections to this process involve parton distributions with four partons. The purpose of this 
paper is to identify and classify these distributions and to calculate their contribution to the 
Drell-Yan process. 

Higher twist corrections, as is well-known, are notoriously complicated. Nonetheless, 
if one is to proceed beyond the naive parton model, it is necessary to identify all parton 
distributions which enter at higher twist, to study their symmetry properties, and to calculate 
their contributions to different hard processes. In contrast to simple distributions involving 
only two quark or gluon fields which occur at leading twist, higher twist distributions involve 
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the matrix elements of the multi quark and gluon field operators at equal light-cone times. A 
complete catalogue of unpolarized and longitudinally polarized distributions up to twist- four, 
and transverse distributions up to twist-three, can be found in ref. 3. 

The results of this paper, which deals exclusively with the twist-four distributions of a 
transversely polarized nucleon, are summarized below: 

Parton distributions are introduced through parton-hadron vertices. As is shown in Fig. 
1, the twist-four distributions in the light-cone gauge are contained in the vertices up to 
four partons. Here the four gluon vertex is absent because it cannot mix with the chiral-odd 
vertices shown. The two-quark vertex, shown in fig. la, provides one twist-four distribution 
h^{x). There are four distributions, named di{x,y) with i = 1,2,3,4 , associated with the 
two-quark-one-gluon vertex in fig. lb. The two-quark-two-gluon vertex, fig. Ic, gives rise to 
three distributions named herein as Di{x, y, z) with z = 1, 2, 3 . Finally, the four-quark vertex 
(fig. Id) has two associated distributions, Wi{x, y, z) with i = 1, 2 . The tensor structure which 
gives rise to each distribution is given in the text, together with restrictions obtained from 
the requirement of PCT invariance. The QCD equations of motion on the light-cone impose 
further restrictions by specifying relations between distributions involving more light-cone 
momentum fractions and distributions involving fewer. Including the distributions identified 
in ref. 3 — i.e. hi {x) , gri^) , Gi {x, y) , and G2 {x, y) — a complete set is now available to deal 
with any hard process involving transversely polarized nucleons up to and including twist- 
four. As one application, we have calculated the twist-four part of the Drell-Yan cross section 
for transversely polarized nucleons. This could provide a framework to analyze corrections if 
an experiment to measure hi{x) is actually performed^. 
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II. TWIST-FOUR DISTRIBUTIONS IN A TRANSVERSELY POLARIZED 
NUCLEON 

In this section, we consider the polarization-dependent twist-four distributions in a trans- 
versely polarized nucleon. To establish notation, we consider a nucleon moving in the z 
direction with its spin vector in the x — y plane with momentum = -p^ -\- ^M'^n^ and 
Sj^ = (0, S±,0) , where S± ■ S± = 1 and p^^ and n^^ are null vectors satisfying = = and 
p-n = 1 . A distribution of partons in such a nucleon is defined by the matrix element of quark 
and gluon field operators at equal light-cone times. We choose the light-cone gauge, ^ • n = 0, 
Simple dimensional reasoning enables determination of the number of inverse powers of the 
hard momentum which accompany a given distribution, and thus its twist. Application of 
PCT rules out a large number of otherwise possible structures. For further details the reader 
is referred to refs. 2 and 3. 

At the level of twist-four, the two quark vertex in fig. la contains only one transverse 
distribution, 

^^^"^^ = ^ / ^'^'^ ^^^^ 1^^°^^^ ^^^""^l ^^"-^ ' 
where A is a soft scale. It can also be projected out from the quark density matrix, 

M{x)= I ^e'^^ {PS^\^{Q)i;{Xn)\PS,_) 

\ (2) 
= -K^i^^$Jiz{x) + ... , 

where dots represent other distributions which are our concern here. 

The two-quark-one-gluon vertex, fig. lb, is more complex and involves two light-cone 
fractions since the quark and gluon can be removed from the hadron from different light-cone 
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"positions" x = Xn or y = jin . There are four PCT allowed transverse distributions: 



di{x,y)= [dXdfx]{PS^\il;{0)iD^{fxn) ■ S^^5ii;{Xn)\PS^) , 



'il;{0)iDj_{ij,n) ■ S±j5-{^ - 7^f)tl;{Xn) 



d2{x,y)^ l[dXdi^](^PS^ 



PS 



dsix^y) = J [dXdfi] {PS^ \iP{0)iD^{fin) ■ iT^iP{Xn)\ PS^) , 
d4{x,y) = j[dXdiJi] /pSi_ 



tl;{0)iDT{fin) ■ iT^-{i4 - i]i:^)ip{Xn) 



PS± 



where — e"^'^^Pi3n^S±s is a vector orthogonal to 5" with T± ■ T± = —1 and, 

The di{x,y) may be readily projected out of the quark-gluon density matrix, 
M''{x,y) =2A2 [ [dXdfji] {PS± \^|;{0)iDl{|ln)^|;{Xn)\PS±) 



Sl^,di{x,y) + Sl-{^ - #)75ci2(x,2/) 



(3) 



(4) 



- ind3{x,y) - iTl-{j>i - ii>)d^{x,y)\ + ... 

It can be easily shown that the diS are real and obey the symmetry relations: 

di{x,y) = -di{y,x) , d2{x,y) = d2{y,x) 



(5) 



(6) 



dz{x,y) = -ds{y,x) , d^{x,y) = d^{y,x) . 
Note that the index a in eq. (5) refers to transverse polarizations (o: = 1,2) only because 

A'^ = and A~ is a dependent quantity whose presence would make a distribution belong to 

one higher twist. 

The two-quark and two-gluon vertex, fig. Ic, also with transverse gluons only, has three 
light-cone fractions. There are three allowed distributions: 

Di{x,y,z) = J[dXdiidiy] {PS± \'{l^{0)iD^{iyn) ■ iD^i/in) i 75 ^{Xn)\ PS±) , 

D2{x,y,z) = / [dXdii du]/pS±\ij{0)[iD±{un) ■ S±i I^^ {iin) + il^^ {un)iD±{iin) ■ S± 



— iD^[un) ■ iD±{iJ,n) ^ j_] ^ 'y5ij;{Xn)\PS± 
D3{x,y,z) = j[dXdiidv\ {PS_l \'{p{0)ie'''^^^p^nsiD_i_a{i^n)iD^p{i^n) i$ ^ P5±> (7) 
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where, 



The above distributions can be projected out of the qGGq density matrix, 



M'^^\x,y,z) ^2K^ j[d\dndu] {PS± \'i/;{0)iD''{un)iD^{^n)ij{Xn)\PSA_) 
+ I {Sll^ + sir - gf ^a) ^ l5D2{x, y, z) 

+ ie'^'^'^^p^ns Ds{x, y, z)j + ... (9) 
Working from the definitions 7a- 7c, certain symmetry relations can be estabhshed, 

Di{x,y,z) = Di{z,y,x), 

D3{x,y,z) ^ -Ds{z,y,x) . (10) 

Finally, the four quark matrix element in fig. Id can be similarly analyzed. The twist-four 
transverse spin distributions are: 

W^ix,y,z) = y y [dXdfidu] {PS^ \iP{0) ^ iP{un)iP{fin) ^ $^ ^{Xn)\PS^) , 
W2ix,y,z) [dXdfidu] {PS^ \iP{0)i^5i^{un)iP{fin) ^ ^{Xn)\PS^) . (11) 

The four-quark density matrix, which is a matrix in two sets of Dirac indices, from which 
these can be projected is, 

M{x,y,z) =2/A^y" [dXdfidv] {PS± \ip{0)^p{un)4;{i^n)^p{Xn)\PSi_) 

=^A^[{^){^, ^J)W,{x,y,z) + {l5m^±)W2{x,y,z)] + ... (12) 
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At this point one needs to take stock of the situation. So far we have, using PCT, Lorentz 
invariance, and dimensional counting, identified all tensor structures needed to describe a 
transversely polarized nucleon at the twist-four level. However, we have further constraints 
available to us in the form of QCD equations of motion on the light-cone, 

V'i = ^T^T^V' • (13) 

This can be used to establish relations between distributions with three light-cone fractions 
with those having two and one, etc. A particularly useful set of relations, which shall be 
needed for establishing the electromagnetic gauge invariance of the Drell-Yan cross is, 

J dzD2{x,y,z) = -y{di{x,y) + d2{x,y) + ds{x,y) + d4{x,y)) , 
J dzD2{z,x,y) = x{di{x,y) + d3{x,y) - d2{x,y) - d^ix^y)) , 
J dydzD2{y,x,z) = x'^h3{x) . (14) 

All distributions discussed so far are implicitly dependent on the renormalization scale. 
Their evolution through radiative processes, however, is a complicated issue and beyond the 
scope of this paper. 

III. THE TRANSVERSELY POLARIZED DRELL-YAN PROCESS 

We now consider application to the Drell-Yan production of lepton pairs in transversely 
polarized nucleon-nucleon collisions. Calculations beyond leading twist require considerable 
formal development, as in the work of Ellis, Furmanski, and Petronzio, and more recently by 
Qiu and Sterman, and Jaffe and Ji.^'^'^ A coUinear expansion of the parton momenta is carried 
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out, different Feynman diagrams are combined together to arrive at colour gauge invariance, 
etc. However, a simpler set of rules emerges from the formalism, which is summarized in 
ref. 3. These rules may be readily applied to the Drell-Yan process by first drawing the set 
of all diagrams which give rise to a lepton pair in the final state, and which do not contain 
more than four partons belonging to a single nucleon. The set of diagrams contributing up to 
twist-four are given in figs. 2-5. All diagrams of a given topology must be included, although 
only one has actually been shown in the figures. 

The hadron tensor that appears in the inclusive Drell-Yan cross section can be written 

as. 



where x{> 0) and y{> 0) are the momentum fractions carried by quarks or antiquarks from 
the nucleon A and B, respectively. For our purpose, we are interested in only the polarization- 
dependent part of the tensor. In the following discussion, we assume that a quark from A 
annihilates with an antiquark from B. for the opposite case, an antiquark from A annihilating 
a quark from B, is obtained by substituting x — > —x and —y — > x. 

To give an example of the application of the rules in ref. 3, we calculate the contribution 
from fig. 2a, 



The factor (—1) arises from the anti-commuting nature of the fermion fields, and an im- 
plicit trace over colour is understood. Ma{x) and MB{—y) are the transverse-polarization- 
dependent part of the quark density matrix. 




(15) 



^''''{x,y)\2a = {-l)TrYMA{x)rMB{-y) • 



(16) 




^75 -^.L hs{x). 



(17) 
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Working out the fermion and colour traces, we find 



+ P^P!ihi{-y)ht{x) 
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Pa-Pb 



S^" — S±A ■ S±B 



^A±B_J^_^_Af_B_ 

Pa-Pb 



where 



S^'^ — Sj^j^Sj^B + ^±A^±B ~ 9t'S±a ■ S±Bi 



(18) 



(19) 



The first term contains only the chiral-odd distributions and the second the chiral-even dis- 
tributions. Quark helicity conservation for massless quarks dictates that the Drell-Yan cross 
section consists only of the combinations of chiral odd - chiral odd distributions or chiral even 
- chiral even distributions. Cp = {N'^ — 1)/2N = 4/3 is a colour factor for three colours. 

Fig. 2b contains one gluon coming from the nucleon A, and the corresponding expression 
for n^"' is, 



YM%{z,x)rMB{-y)ilo 



(20) 



[z - x) i>j^ -y i>B_ 

where the twist-four part of the matrix Ma{z^ x) is given in eq. (5), and the twist-three part 
in eq. (43) in ref. 3. After a lengthy calculation, we find the contributions from all four 
diagrams having the fig. 2b topology, 

■''"<-«Ww + ^.lf(-,)l 



y 



X 



y 

+ / dw —Dii-y, -z, -w) 

X I z 



dw — D2{x^ z, w) 



^ ^A^A + y ^B^B 



— -jCpA^iS^LA ■ S±b) 
+ IcfA'{S±a ■ 8^3)9^" 9?i-y) I dz 



9T{-y)9T{x) 



xyPA ■ Pb 

X — z 



9t 



dz 



Gn-y,-z)-G§{-y,-z) 
y- z 



(21) 



The first two terms in the above formula are chiral odd, and the second two are chiral even. 
The transverse metric tensor has only = = — 1 as non-zero components. 

For the diagrams with two gluons from the same nucleon shown in figs. 3a and 3b, there 
are only chiral-odd contributions. Using the density matrix in eq. (9), we find. 
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^^''{x,y)\sa^lCFA\S^A-S^B)[x'P!^Plhf{-y)hi{x) 



1 



n^''{x,y)Ub = --CFA'S^'' 



+ y''P^PBhi{-y)ht{x) 

hf{-y) 



1 



(22) 



x^y^PA-Ps' 

D^{x, X2, Xi) + D^{x, X2, Xi] 



+ ^ ^ ' I dxi dx2 



dxi dx2 

y J x-xi 
Di{-y, -XI, -xx) + Df{-y, -X2, -xx) 



(23) 

X J y — xx J ' 

In eq. (22), we have made use of eq. (14c). Clearly, the chiral-odd longitudinal part from eq. 

(22), combined with these from eqs. (18) and (21) is electromagnetically gauge invariant. 

There are four diagrams in which one gluon comes from A and the other comes from 

B. These diagrams contribute only to chiral-even part of the cross section. Our calculation 

shows. 



fl^^''{x,y)\u =CfA\S^a ■ 8^3)9^" j ^^^^ [Gi{x,w) - G^{x,w)\ 
X [Gf (-^, -y) - Gf{-z, ~y)] +{x^-y,A^ B) 



(24) 



1 



dz dw 



--CaA'{S±a-S^b)9' 



Gt{x,w)-Gt{x,w)\ 



X 



Gf{-z, -y) - Gfi-z, -y) \^{x^-y,A^B) 



(25) 



where Ca = 3. The second term in eq. (25) requires some explanation: Gi{x,w) is defined 
exactly as Gi{x,w) as in eq. 38 of ref. 3, except that the covariant derivative has only the 
gluon potential. This is because the partial derivative coming from transverse momentum 
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expansion of fig. 2a does not contribute to terms multipfied by Ca- Since the final result 
has to be gauge invariant, we replace the by in~ F^-^ / {w — x), and thus, Gi{x^w) is 
manifestly gauge invariant. And, 



J z w 

X [Gf{-y,-z)+Gi{-y,-z)] 



dz dw 



2 J z w 



G^{x,w) + G2ix,w) 



X 



Gf{-y,-z) + Gn-y, -z) 

dz dw 



(26) 



n^''{x,y)Ud = -CAA\S^A-S^B)g^'' I -^[Gt{x,w)Gf{-y,-z) 

J z — y w — X 



G^{x,w)G^{-y,-z)\ 



(27) 



Again, it is simple to see that the chiral-even longitudinal part in eq. (24), combining with 
these from eqs. (18) and (21) is electromagnetically gauge invariant. 

Finally, we consider the four-quark diagram shown in fig. 5. Evaluation requires a double 
trace as there are two quark loops. This is straightforwardly done, and using eq. (12), we get: 



n^''{x,y)\^ = 



\c^^2gt,v^}i{-y)_ f dz dw 
4 2y J z — X w — X 



J2[w^'{x,z,w) 



+ W^^{x - z, -z, w-z) - Wf{x, z,z-w) - Wf'ix - z, -z, -tu)] 



l^^,^^2^M^ ^i (^ ) f dz dw 



z — y w — y ^ 



+ W,^{x- z, -z, w-z)- W,^{-y, -z, -w) 



W,^{-y + z,z,-w + z) 



(28) 



The different terms in eq. (28) comes from eight different diagrams having the same topology 
as fig. 5. 
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It may be readily verified that our final result for W^^'' , which is given by the sum of all 
contributions from diagrams 2-5, is electromagnetically gauge invariant: q^W^jy = W^^^q^ = 
. This is an important partial check on the correctness of the calculation. 

IV. SUMMARY 

In this paper, we introduce the polarization-dependent twist-four distributions in a trans- 
versely polarized nucleon. The most general ones contain three light-cone momentum fractions 
and are Di and Wi defined in eqs. (7) and (11). There distributions, together with the twist- 
two and twist-three distributions defined for the same nucleon state, form a complete set for 
describing any hard scattering processes involving traverse polarization. 

As an example, we have calculated the twist-four correction to the Drell-Yan process. 
The final result is gauge invariant in both electromagnetic and colour interactions. For lack of 
information on these distributions, we cannot access numerical significance of the correction 
other than simple dimensional analysis. However, our formulas can be coupled with model 
calculations of the distributions to give a detailed prediction. 
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FIGURE CAPTIONS 

For a hard copy of the figures, please send email to ereidell@marie.mit.edu 

Fig. 1: Parton density matrices which can contribute to a hard process at the twist- four level. 

Fig. 2: A quark from hadron A annihilates with an antiquark from hadron B leading to a large 
mass virtual photon. Other diagrams of this type are obtained by allowing the gluons in 
(b) to interact with the quark line on the right, and interchanging A and B. 

Fig. 3: Representative diagrams which give chiral odd contributions to the twist-four Drell-Yan 
cross section. 

Fig. 4: Diagrams involving one gluon from each nucleon which give chiral even contribution to 
the Drell-Yan cross section. 

Fig. 5: Four quark diagram contributing to the Drell-Yan twist-four cross section. 
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